Predicted defect induced vortex core switching in thin magnetic nanodisks 
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We investigate the influence of artificial defects (small holes) inserted into magnetic nanodisks 
on the vortex core dynamics. One and two holes (antidots) are considered. In general, the core 
falls into the hole but, in particular, we would like to remark an interesting phenomenon not yet 
observed, which is the vortex core switching induced by the vortex-hole interactions. It occurs for 
the case with only one hole and for very special conditions involving the hole size and position as 
well as the disk size. Any small deformation in the disk geometry such as the presence of a second 
antidot changes completely the vortex dynamics and the vortex core eventually falls into one of 
the defects. After trapped, the vortex center still oscillates with a very high frequency and small 
amplitude around the defect center. 

PACS numbers: 75.75.+a, 75.60.Jk, 75.60.Ch, 75.40.Gb 



I. INTRODUCTION AND MOTIVATION 

As it is well known, magnetic vortex states are exper- 
imentally observed in ferromagnetic disk-shaped nanos- 
tructures. These topological objects exhibit a planar- 
like configuration of spins outside the core, where a per- 
pendiculajpiptnagnetization (up or down polarization) is 
observecocrB. As long as one could manipulate these 
vortices, several possibilities would emerge. Therefore, 
studies taking into account the application of external po- 
tentials such as static and sinusoidal ma^a^C|tic Jelds are 
nowadays very common in the literatureQ'ElaEHa'EI. For in- 
stance, the simplest effect induced by an external field is 
the gyrotropic mode, which is the lowest excitation of the 
vortex structure. This mode is simply the elliptical vor- 
tex core motion (with the resonance frequency) around 
the disk center. The sense of gyration in an elliptic 
trajectory (clockwise or counterclockwise) is determined 
only by the vortex core polarization. External potentials 
such as magnetic fields or (d.c.) spin-polarized currents 
can also S|tiioulate dramatic effects such as the switching 
behaviorElBBEl. However, to the authors knowledge, no 
internal mechanisms able to induce the core switching 
were reported neither experimentally nor theoretically. 
Such a chance would only occur if the vortex could inter- 
act with possible inhomogeneities present in the nanos- 
tructure. For instance, a possibility of triggering this pro- 
cess could be obtained by removing some small portions 
of the magnetic nanodisk, in such a way that the cavities 
(antidots) so created work ^^^|Pfi^\VSi"Ii"^i eventually 
affecting the vortex s^uctu Recently, 
Rahm and coworkersli30li30 have experimentally stud- 
ied the cases of one, two, three and four antidots arti- 
ficially inserted in a disk with diameter ^ 500 nm and 
separated by around 150 nm — 200 nm. Their results 



not only confirmed the previous stateraejit about vortex 
pinning around the hole defectslijoli3'E3 as well as have 
put forward the possibility of using these stable states 
as serious caiididates for magnetic memory and logical 
applicationgla. Although experimental results are pro- 
vided for such systems, it is worthy noticing that a broad 
theoretical analysis is still lacking. In this paper, we ar- 
gue that defects present in the dots may induce interest- 
ing possibilities for the magnetization dynamics, includ- 
ing the switching behavior of vortex structures. This is a 
very interesting topic and open new possibilities for teclt 
nologies and experiments. Indeed, a very recent worklij 
has shown that magnetization dynamics and evolution is 
deeply modified when cavities are introduced. Their sam- 
ples consist of Co-made elongated ring-shaped nanomag- 
nets, which like Permalloy has negligible magnetostatic 
anisotropy. What is clearly observed in these experi- 
ments is that the reversal pathway of the so-called diago- 
nal onion state is drastically affected by the inserted holes 
once they attract and sometimes pin the transient vortex- 
like domain- wall configurations. Actually, not only the 
presence, but also the sizes and locations of the cavi- 
ties (two or four in their hollowed samples) remarkably 
change magnetization evolution and eventually lead to 
its reversal. We will follow a different line in which the 
magnetization reversal takes place due to the effects pro- 
duced by the hole on the vortex core motion. No further 
agent is necessary in the switching process, like occur to 
other usual methods where external fields are one of the 
main ingredients. 

Before starting our analysis, it would be useful to 
describe how a magnetic vortex can arise in nanodisks. 
Indeed, a small revision will help to justify the magnetic 
model used here. The vortex state in magnetic nanodot 
materials is the result of the competition between ex- 



change and magnetostatic interactions. Particularly, : 
finite systems like magnetic dots, densities of "magnet 
charges" are induced in their volumes and surface 
Considering the magnetization M, the volumetric ar 
superficial densities are defined as p = — V • M ar 
a = M ■ h respectively, where n is the unit vecto 
normal to the surfaces of the dot at each point. Ther^ 
fore, in order to minimize the magnetostatic cncrg 
the magnetization tends to point parallel to the d< 
surfaces and it is one of the properties responsible fi 
the formation of the vortex ground state. Howeve 
as long as we go towards the vortex center, exchanj 
energy density gradually increases like r~^. Then, : 
order to regularize exchange contribution, the magnet 
moments around the center tend to revolve, dcvclopir 
out-of-plane projection, so that exactly at the vort( 
center it is perpendicular to the disk face, whoi 
direction defines what is called vortex polarizatic 
(p = +1, if the centered moment points up, or p — — 
if it is pointing down). Here, we study very th: 
ferromagnetic nanodisks with thickness L and radii 
R so that the aspect ratio L/R <c 1. In this case, oi 
can assume V • M = and hence, the only source i 
magnetostatic interactions is the superficial magnetic 
charges. To write the vortex state in these dots, it is 
convenient to parametrize the magnetization Af(r) by 
two scalar fields, the polar 9(f) and azimuthal (/){f) 
angles, j2{r) = M{f)/M = (sin cos (/i, sin sin 0, cos 0). 
Using m(r) = cos^?(r), the vortex state configuration 
can be written as m(r) = 'rn^i'r) = cos 9q {r) , 4>{r) = 
(r) — arctan(2//x) ± tt/2, where the magnetization 
unit vector deviates out of the plane at the vortex core 
™o (0) ±1 and outside this region, jtiq (r) — > 0. 
The vortex core size is approximately the exchange 
length Zo (which can be taken to be of the order of 
the lattice spacing a in our simulations on a discrete 
lattice). In this notation, the values of m^(0) at the 
vortex center defines the vortex polarization. On the 
other hand, (^^(r) defines the vortex chirality. The flux 
closed in-plane vortex {(j)Q{f)) and the perpendicular 
component in the center (m^(())) are independent of 
each other so that four different magnetization states of 
one disk are possible. 



Our paper is organized as follows: in Section II, the 
model and associated methods of numerical analysis are 
presented and discussed. Section III is devoted to the 
study of sample with a single hole. There, we emphasize 
the vortex core reversal process induced by its interac- 
tion with the defect and the non-ordinary circumstances 
favoring it. The case of doubly hollowed disk is studied 
in Section IV, where differently from the former case, no 
holes configuration yields the core switching: indeed, it 
is always captured by one of the defects before any re- 
versal be observed in our simulations. We finally close 
our work by point out our Conclusions and Prospects for 
forthcoming investigation. 
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FIG. 1: Trajectory of the vortex core (beginning at the disk 
center (20a, 20a)), in the xy-plaxie, after turning off the ex- 
ternal field. Here, no hole is inserted and the arrows indicate 
the direction of the motion. 



II. THE MODEL AND METHODS 

In order to describe a thin magnetic nanodisk we pur- 
sue an alternative strategy substituting the disk by a two- 
dimensional film with magnetization vectors distributed 
in a regular square lattice inside a circumference of radius 
R. In addition, the magnetostatic interactions due to the 
presence of the magnetic charges in the lateral and top 
surfaces of the dot are replaced by local potentials. Since 
we are going to study systems with one and two holes 
and also with no hole, the model used is summarized by 
the following general Hamiltonian: 

-H" = - X] -^iil^i ■ '"j + ^^aiflk ■ nk,a)'^ 

a=l,2 k 

+ X! '^^DAfll-f>-D,lf-^h-ilt, (1) 
D=l,2 ; i 

where Jij = for sites inside the holes and Jjj = 
J > for the remaining sites of the film. Here, fli = 
Mi{r)/Ms = i^ifx + fi\y + jifz is the unit spin vector 
at position i (Ms is the saturation magnetization), the 
sum {i,j} is over nearest-neighbor spins, the terms with 
positive constants {Xa,^D,i) mimic the magnetostatic en- 
ergies (see below) in the top face of the disk, Ai, while 
A2 accounts for the lateral edge. At the holes borders 
the analogous constants taking into account the surface 



charges at the edges are Xi,i and A2,i, for holes 1 and 2 
respectively. For the case without defects, Ai^; — \2,i = 
while for the situation with a single defect A2.; = 0. In 
turn, h is an external magnetic field. The hole defects 
are introduced by removing a number of neighbor spins 
around sites fi = (0,di) and r2 = (0,-^2) from the sys- 
tem (di > and c?2 > 0). Of course, the number of 
neighbor spins removed around a particular position (ri 
or r2) defines the hole sizes qi and Q2- Now, the mag- 
nctostatic energy favors the spins to be parallel to the 
film surfaces. Therefore, the local unit vectors h defined 
on the points characterizing the surfaces of the material 
are always perpendicular to these surfaces (or bound- 
aries): nfc 1 = z IS perpendicular to the disk face, in the 
xy-plane; fik^2 is radially perpendicular to the circumfer- 
ence envelop, with radius R, at each point; fiij and 7^2,; 
are radially perpendicular to the border points of holes 1 
(pointing to fi) and 2 (directed to r2), respectively. 

As we have already mentioned, in principle, our model, 
Hamiltonian (|]), could describe a very thin disk with 
aspect ratio L/R <C 1. In fact, several known experi- 
mental results with nanodisks can be qualitatively repro- 
duced using it. For instance, the vortex core pjcecession 
mode (which has a rather low eigenfrequencyocj) can 
be easily obtained, at least qualitatively, with this sim- 
ple model (see below). However, it is not well known 
how defects can change this picture. De£c|jts_Hiust be 
naturally presenttil or artificially insertedll300_ij|L^h£i 
film samples. It has been predicted.|theoreticallyll3lijll3'E^ 
and also observed experimentallyli3n3 that small holes in- 
corporated into the magnetic structures can attract and 
capture magnetic vortices. On the other hand, we argue 
that another effect can occur: the vortex core switch- 
ing due to a involved dynamical state of the magnet, 
which should include the coherent magnetization oscil- 
lations (near the defect) matched to the vortex motion. 
Actually, a somewhat similar effect has been observed 
in quite recent experimentsEEI. There, the pathway de- 
veloped to reverse the diagonal onion state in elongated 
nanorings is quite sensitive to intentionally inserted cav- 
ities, as well as their sizes and locations. Here, the pic- 
ture comprises a vortex-like configuration with polariza- 
tion, say, up, which after interacting with a hole, of a 
given size at a specified position (these values clearly de- 
pend on other parameters like the whole nanodisk size, 
applied field amplitude for resonantly exciting the gy- 
rotropic mode, etc), turns it to be down, eventually re- 
versing the sense of its gyrotropic motion. What should 
be stressed is that such a remarkable phenomenon show 
up in our simulations only for very special situations in- 
volving the physical parameters, for instance, the hole 
must be relatively small (around 5%-10% of the sample 
size). In addition, the insertion of a second similar hole 
in any position deeply jeopardize the magnetization evo- 
lution and the vortex core is always captured by one of 
the holes before any core reversal process. 

Our results are obtained by using spin dynamics sim- 
ulations for a square lattice occupying all possible points 
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FIG. 2: Magnetization in the a;-direction as a function of 
time and its Fourier transform. 



inside a circumference of radius R (in most simula- 
tions we have used R = 20a, 25a, 30a). It is solved by 
employing the fourth-order predictor-corrector method. 
The calculations presented here consider only R = 20a, 
Ai = 0.2J and A2 = 2J. The choice of other values 
of Ai does not alter the essential physics if it is not 
large enough (in all investigations we have used the range 
< Ai < 0.28). For Ai > 0.28, the vortex becomes es- 
sentially planar and does not develop the out-of-plane 
components at the corecil. On the other hand, we have 
verified only a limited range of A2. Indeed, the inter- 
val 1.8 < A2 < 2.2 leads to the same quantitative and 
qualitative results. In order to excite the vortex core 
precession mode, an external perturbation was applied. 
Here, we have chosen a sinusoidal external magnetic field 
of the form hext{t) = /iosin(a;t), with = O.OlJx and 
bj = 0.089 J. These are the best conditions we have iden- 
tified in our simulations favoring the core reversal. For 
instance, large changes in these field parameters will not 
produce the main phenomena we intend to describe. In 
all cases, the duration of the sinusoidal field was 700 
and the total time of the simulations was 10^J~^ with 
time increment equal to At = O.OIJ^^. Firstly, in the 
absence of the defects and after turning the field off, we 
observe the precession mode with a circular trajectory, 
as expected (see. Fig. ^. For a better analysis of this 
mode, we have plotted the behavior of the magnetization 
in the x-direction < ^l^ > as a, function of time and its 
Fourier transform (see Fig. ^; the same behavior can 
be seen for < /Zj, >). The Fourier transform presents a 
peak at a well defined frequency ojn = 0.0056J, which is 
the frequency of the vortex core precession (the resonance 
frequency, see Fig. ||). Of course, it depends on the size R 
of the disk. More realistically, the resonance frequency 
depends on the aspect ratio L/R. However, since our 
model assumes a very thin disk, we can, in principle, ob- 
tain LOn as a function of L/R only ior L/R << 1 . We also 
studied the behavior of ujji as a function of 1/i?, which 
appears to be essentially linear, at least for 1/R << 1 
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FIG. 3: Position components {cx,Cy) of the vortex center as 
function of time. Note that the Fourier transform leads to a 
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FIG. 4: The behavior of the gyrotropic frequency, ujr, as 
function of the 'aspect ratio', R~^, obtained from the present 
model. Points are simulational results while the line is the 
linear fitting. 



(see Fig. HL in. qualitatively accordance with experimen- 
tal findinga3'EHl. All above comparisons (and others not 
presented here) show that this model works very well. It 
gives us a perspective that Hamiltonian (Q) is also able 
to predict new facts. 

Before studying the problem of interest, we have also 
plotted the behavior of the z-component of the magne- 
tization as function of time (and its Fourier transform) 
in the absence of a defect (see Fig. |^). In this case, the 
average magnetization over the whole sample oscillates 
around a small positive value (because the out-of-plane 
spins forming the vortex core are pointing along the pos- 
itive z-direction) all the time without abrupt changes, 
indicating that there is not any magnetization reversal 
during the vortex core precession, as expected. In addi- 




FIG. 5; < j-iz > versus t and its Fourier Transform. Note 
that < > oscillates around a positive value. 



tion, the main peak in the Fourier transform occurs at a 
frequency 0.078J, which is different from ujj^. 

This model involves several parameters and, of course, 
it is very hard to treat in details all aspects. Therefore 
we have exposed here the general behavior of the vortex 
core dynamics and the main particular phenomena 
found for specific conditions. We notice that a hole, 
independently of its size and position, in general, at- 
tracts the vortex core to its center, i.e., the core almost 
always falls into the hole (sometimes , | pcjat j U'r |i ng -like 
events may occur). Some theoreticalll3llJE3'E3E3 and 
experir HCQt al . works have already reported the capture 
processEllESEZl. Hence, it would be interesting to consider 
other possibilities for the vortex dynamics in the pres- 
ence of such type of defects. Two situations are analyzed 
in details. The first one, which is more remarkable, 
describes the vortex core switching due to the vortex- 
defect interaction. It happens only for very special 
circumstances and our work indicates the most suitable 
configurations, in order to aid possible experimental 
probing. The another one relies on the fate of the core to 
be trapped by one of the two holes (two configurations 
that yield this capture are presented). Such a situation 
also serves to illustrate how a relatively small pertur- 
bation eliminates the magnetization reversal and how 
the center of the vortex oscillates after a capture process. 



III. A NANODISK WITH A SINGLE HOLE 
AND THE VORTEX CORE SWITCHING 

Now we consider some effects that a hole may induce on 
the vortex core dynamics. When a small hole is inserted 
in the circular film, the vortex core suffers new influences 
so that the magnetization dynamics should become much 
richer. Indeed, ho lc a i ftduce an attractive effective poten- 
tial for the vorte3al3llil'E3'Ea and therefore, such force may 
compete with the effects of applied fields, magnetostatic. 
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FIG. 8: Magnetization in the «-direction and its Fourier 
transform. The vortex-hole interaction causes the switching. 
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FIG. 6: Typical trajectory of the vortex core in a disk con- 
taining a hole (at position (0,33a)). The vortex core motion 
starts in the disk center (20a, 20a) and the precession mode is 
excited by a field pulse. The comparison of this motion with 
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FIG. 7: < fix > versus t and its Fourier transform for a film 
with a hole. 



anisotropies, etc. 

In this Section we consider the presence of a sin- 
gle hole of radius gi ~ g. For the external field, we 
use the parameters of the previous Section. In a nan- 
odisk with a cavity, one expects that magnetic charges 
are also induced in the cavity walls analogously as they 
appear in the usual surfaces of the disk. Therefore, tak- 
ing into account only the presence of defect 1 at posi- 
tion ri, in principle we should use Ai^; ^ 0, which cor- 
responds to the influences of the magnetostatic energy 



coming from the edge of the hole cavity. Here we will 
show only influences of very small defects as compared 
with the disk size and consequently, Ai^; could be ne- 
glected in a first approximation since the hole wall can 
support only a tiny amount of magnetic charges. More 
realistic analysis should then take this parameter to be 
small, for example, as compared with A2. However, it 
is illustrative to study the system for a large range of 
Xi,i values, say, < Ai_; < A2, for elucidating its effects 
on the whole properties of the magnetization. Gener- 
ally, for the smallest hole considered (four neighbor spins 
removed, gi ~ g ~ a), the vortex core is almost al- 
ways captured (or sometimes scattered) by the defect for 
Xi.i > 2.3 X 10-4. Hence, large enough contributions of 
magnetostatic energy coming from the cavity wall (i.e., 
large Xi,i) leads to the common phenomena of capture, or 
even the more seldom scattering-like effect, observed for 
larger holes. Indeed, larger holes always capture the vor- 
tex core, independently of Ai_; (recall that besides mag- 
netostatic, exchange gets loweij-as-tlie hole size increa se s ) 
what agrees with experimentaltilOEO and theoreticalB£lE3 
results, which consider defects occupying an appreciable 
fraction of the disk. Therefore, if one wishes looking for 
different possibilities, perhaps small defects should be the 
focus. With this strategy, we first analyse several pa- 
rameter possibilities and later we shall consider g — a 
and use Ai,; — A2/IO4 (or, for simplicity, Ai,; = 0) for 
all subsequent calculations. Indeed, the most important 
phenomenon we would like to describe here happens only 
for < Ai.i < 2.3 X IQ-^J. Physically, it is reasonable to 
think that a small cavity wall has only a tiny contribution 
to the magnetostatic energy. 

The size and position of the hole can be controlled 
and so we have studied several cases. For example, for a 
central hole with size a, the resonance frequency of the 
gyrotropic mode is given by 0.0034J, which is smaller 
than the one we have observed for the usual nanodisk 
LOn = 0.0056<i-(a fact qualitatively predicted by analyti- 
cal treatmentH). The problem with this configuration is 
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that the ground state is a vortex pinned to the hole (of 
course, without an out-of-plane core) and, therefore, the 
vortex needs to be previously released from the cavity 
to exhibit the gyrotropic mode. Changes in the defect 
size and position lead to the following general situations: 
depending on the initial conditions (external sinusoidal 
fields, essentially), then, as the vortex core moves suf- 
ficiently close the defect border, it can be either trans- 
mitted or captured by the cavity. These are interesting 
effects due to core- hole interaction but, for very special 
features of the film, a much more rare dynamical event 
may be triggered. Actually, choosing the coordinate sys- 
tem in which the disk center is put at (20a, 20a) and if 
and only if we use fi = (0, 33a) (or (0, 7a) ) for the cen- 
ter position (along the y-axis at di = ±13a) of a hole 
with size g = la (four neighbor spins removed), then a 
fabulous phenomenon occurs: the magnetic vortex core 
is reversed. For all disk sizes studied, only one configu- 
ration was found to be able of reversing the core magne- 
tization. Always, the hole must be relatively small and 
placed near enough the disk border, perpendicularly to 
the field. Particularly, for R = 20a and p = a, the asym- 
metric conditions with di = 13a seems to be the unique 
defect position possible to trigger the switching process. 
It should be emphasized that the switching of magnetic 
domains generally depends on a myriad of detailed fea- 
tures of the magnetic particles, and topological effects 
ultimately limit this possibility. Therefore, the occur- 
rence of such a delicate dynamical effect should not be 
expected for any ordinary configuration of the magnet 
with a hole. Usually, what is observed is the capture of 
the core when the defect is large enough; conversely, only 
small changes in its dynamics are observed for very small 
defects. The main lesson is: the core reversal seems to de- 
mand relatively small hole (in our case around 5% — 10% 
of the nanodisk dimension) combined with other favor- 
able events, like its location and suitable applied oscil- 
lating field. Another important point is that, assuming 
Ai,i/A2 << 1 and for special conditions and parameters 
(like di and g used here), the reversal of the magneti- 
zation of the vortex core (along the z-direction) can be 
always reproduced. In Fig. ^ we show the trajectory of 
the vortex core. Note that this mode starts with the 
same frequency of the case without a hole (Fig. |l|) in a 
trajectory outer the hole. However, after some revolu- 
tions, when the vortex core moves again very near the 
hole border (in a region in between the hole and disk 
walls), it quickly changes the direction of the motion si- 
multaneously with the magnetization reversal along the 
z-direction (see video of the vortex core motian and the 
switching phenomenon as auxiliary material23). Really, 
a change in the sense of gyration unambiguously indi- 
cates a change in the vortex core polarization. Then the 
vortex core again develops a new precession mode with a 
smaller frequency = 0.0025 J as shown in Fig. |^ (an 
orbit with radius smaller than dy). An appreciable part 
of the kinetic energy of the vortex was transformed in 
spin waves. To see more clearly the switching behavior, 




FIG. 9: (Color online) The instant immediately before the 
switching process. The out-of-plane magnetization fluctua- 
tions can be perceived in the irregularities present around 
the green surface. Note the deep protuberance with negative 
values of fiz near and below the hole and the core. Such a pro- 
tuberance will pull the vortex core, leading to the inversion 
of its polarity. 



the z-component of the average magnetization < p.^ > 
is plotted in Fig. ^, which can be compared with Fig. 
||. Note the rapid change of the mean value oi < fiz > 
from positive (when the core direction was up) to nega- 
tive (when the core is down). In the Fourier transform, 
the main peak happens in almost the same frequency for 
the case without a hole. 

Thus, the lacking of experimental observation of core 
switching induced by vortex-hole interaction may be 
credited to the large holes generally inserted into the 
samples, around 15% — 25% of its size. In these cases the 
commpn phenomenon is the capture of the core by the 
defectE3E3EZI) . Basic physics concepts may help to under- 
stand this even better: as it is known the effective poten- 
tial (exchange -I- magnetostatic) experie nced , byj-the vor- 
tex is globally minimized inside the holdlSllil'Ea'ES. Con- 
versely, if the hole is very small, it only slightly affects 
vortex dynamics and no appreciable effects are expected 
to occur. Actually, the phenomenon of core reversal pre- 
dicted here takes place for a hole of intermediary size, sit- 
uation in which the minima provided by the exchange and 
magnetostatic contributions are comparable. In sum- 
mary, our simulations as a whole agree with experimental 
observations upon vortex core dynamics in the presence 
of intentionally inserted cavities: generally, the core is 
captured by one of the defects, mainly when the y ar e 
sufficiently large, like in the experiments of Refs.t3'E3EZl. 
Only for special circumstances the core reversal is ex- 
pected to occur, demanding among other a non-centered 
relatively small hole. Another peculiarity of core rever- 
sal induced by vortex-hole interaction is the fact that the 
polarization switching occur only from internal dynam- 
ical events: the vortex must be only under motion (the 
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gyrotropic is the more appropriate for it is easily repro- 
ducible) and pass suitably close to the defect, so that the 
attraction from the hole ultimately dominates its inertia, 
reversing its motion sense. This in turn, demands the 
switching in its polarization, so that the sign of the gy- 
rotropic factor in the Landau-Lifshitz-Gilbert equation, 
G X V is kept unaltered. Another dynamical explana- 
tion for this phenomenon is: firstly, the vortex motion 
excites spin-wave modes. Differently from the case with- 
out defects, the presence of the hole makes these modes 
with larger amplitudes due to reflections in its border, 
mainly when the vortex core moves near it. Secondly, for 
a magnet of finite size, one expects that the radiation of 
spin-waves, their reflection off the disk (and hole) bound- 
aries and their effect back on the vortex core will result 
in the establishment of a dynamical state of the magnet 
which includes both the confined moving vortex and the 
coherent magnetization oscillations matched to the vor- 
tex motion. Then, under this special situation, strong 
out-of-plane spin fluctuations arise near the hole border, 
changing coherently the direction of the vortex core mag- 
netization when it passes close to the defect. Figure ^ 
shows an instant of the core motion immediately before 
the magnetization reversal. 




FIG. 10: An instantaneous position of the vortex core in a 
circular magnetic thin film with two holes centered at (0, 13a) 
and (0, —8a). The holes sizes are of the order of the lattice 
spacing a (in each defect, there are four spins removed). The 
length of the arrows are proportional to the spin projection 
into the XY-plane. Note the small arrows around the vortex 
center at (8a, —2a) forming the vortex core. 



IV. A NANODISK WITH TWO HOLES 



Our aim in this Section is to consider the effect of two 
holes on the vortex core dynamics in circular thin films. 
Also, this study will be useful for showing that the rever- 
sal magnetization phenomenon is completely suppressed 
if a second defect, even small and distant from the first, 
is introduced into the system. In addition, the pinned 
vortex dynamical behavior is detailed. Like the former 
case, the disk radius is selected to be i? = 20a. 

The results are obtained by choosing an external si- 
nusoidal field to excite the gyrotropic mode. Again, the 
field is applied for a short time (~ 700J~^) and its pa- 
rameters are ho = O.OlJi and uj = 0.0093 J. Of course, 
the vortex motion will depend on the holes size and po- 
sition as well as the external field. In all observed cases, 
the vortex finishes trapped in one of the antidots just it 
finds the defect in its way, even for smaller holes. Since 
this is the common phenomenon for all possible positions 
of the defects, we intend to know the dynamical process 
after the trapping. To illustrate this, we choose the fol- 
lowing geometry: hole 1 at position fi — (0, 13a) and 
hole 2 at r2 = (0, —8a) as shown in Figs. ( [l0|) and (11). 
Larger defects will capture the core even faster. What 
should be stressed is that the presence of the second de- 
fect, even if smaller and distant of the first, will change 
the vortex dynamics drastically so that no core switching 
have been observed before its capture by one of the holes. 
Figure ( |l^ shows the mean magnetization < /i^ > along 
the x-direction as a function of time. As the vortex core 
moves, this mean magnetization oscillates almost har- 
monically from —0.5 to 0.5, indicating that the motion 
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FIG. 11; (Color online) Three-dimensional view of Fig. (^o|). 
The vortex polarization is up (yellow-red peak). Note the 
out-of-plane fluctuations generated by the core motion. These 
spin waves propagate through the system and are reflected in 
the disk and holes borders. 



is approximately circular (or elliptical) around the disc 
center. The orbit is outer the holes. However, at about 
t ~ 5000 J~^, the core decreases considerably its velocity 
near the hole 1, moves slowly in spiral (outer the anti- 
dot 2 and inner antidot 1), until hitting antidot 2 again, 
where it is captured at about t = 8000A~^. Now, the 
sense of gyration is not changed and consequently, there 
is no switching process. The < jj,^ > magnetization oscil- 
lations stop abruptly becoming almost constant around 
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FIG. 12: Mean magnetization < /i^ > as a function of time 
and its Fourier transform. The resonance frequency is approx- 
imately 0.0057J. This figure shows clearly that the vortex 
cor< 
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FIG. 13: Mean magnetization < /i^ > as a function of time 
and its Fourier transform. 

0.4 (actually, < fi^ > oscillates very rapidly with a very 
small amplitude around 0.4). We also have calculated 
the Fourier transform of < /i^ > (see Fig.(p^)). It has 
two main peaks: the first one refers to the oscillating 
external field (along the x-direction) and the second is 
related to the resonance frequency. In our unit system, 
such resonance frequency is given by ujres-d ^ 0.0057J, 
which can be compared with the analogous problem of a 
film without defects lor ~ 0.0056 J. Similar graphics are 
obtained for < > as shown in Fig.(p^). The difference 
is that, after being pinned, the mean y-niagnetization os- 
cillates around zero (with a relatively large amplitude as 
compared to its counterpart in x) and not around a finite 
value as happens to < /i^ >. In addition, there is no the 
first large peak exhibited in the < /i^ >. The reasons 
are simple: first, the holes were displayed in the y-axis 
(at (0, 13a) and (0, —8a)) and second, the external sinu- 
soidal field (responsible for the first peak in the Fourier 
transform of < fj,^ >) is applied along the x-direction and 
hence, it can not appear in < /i^ >. The trajectory of 
the vortex core is shown in Fig.(Q). 

Based on the above description, an interesting fact to 
report is that the vortex center oscillates around the hole 
center with small amplitude and large frequency when 



the vortex is pinned to the defect. This phenomenon 
looks like the one thepijetically predicted in layered 
ferromagnetic materialsE3. Nevertheless, this effect 
must be very richer in nanodisks because the surface 
magnetostatic energy has an effective contribution, 
which is— absent in the layered systems. Analytical 
estimates predicts that such oscillations are related to 
the vortex mass and take place around 10^ GHz, which 
is far bejipnd current possibilities of observation in these 
systemsE^I, ~ 10 GHz. However, a study considering all 
microscopic details of the pinned vortex oscillation will 
require further investigation. 

In our present analysis, the oscillation frequency of this 
mode can be seen in the Fourier transform of < /i^ > or 
< > and is given by ~ 0.123 A. It is much higher 
than the resonance frequency, and much lower than mass- 
like vibrations, as above. Such an intermediary frequency 
could be somewhat related to the fractional gyrotropic 
mode, which is predicted to occurs around IO^uir and 
is provided by a remanent (fractional) gyrotropic vea- 
tor of the vortex core, partially captured by the holeE3. 
Another interesting feature can be easily seen when one 
compares the Fourier transforms of < /i^ > and < /z^ >. 
Note that, after being captured, the magnetization am- 
plitude oscillations are larger in the y-direction than in 
the x-direction. It means that the pinned vortex cen- 
ter oscillates much more along the x-axis than along the 
j/-axis. Therefore, it describes an ellipse with larger semi- 
axis along the x-direction, of the order of the hole size. 
For a trapped vortex, the difference in the amplitudes of 
oscillation of < /Lt^ > and < /x^ > increases as the hole 
responsible for the capture is placed nearer the circum- 
ference border of the film. Hence, this ellipse turns out to 
be an almost straight line if the hole is located very near 
the disk boundary, while it takes a circular shape if the 
cavity is centered in the disk. Therefre, several proce- 
dures of production and propagation of spin waves with 
different phases could be built only by changing the de- 
fect position, which is easily manipulated by lithographic 
techniques. Thus, our results indicate that the vortex- 
hole system could also works as a mechanism to create 
and control spin waves in confined magnetic thin films. 

We have also plotted the mean magnetization along 
the ^-direction, < /i^ >, and its Fourier transform. 
Fig. ([Tsl) . There, we clearly realize that, before the 
capture process, < fi^ > oscillates around a small 
positive value, once the vortex is up-polarized (negative 
if the polarity were down). After some vortex core 
laps, at about t ^ SOOOJ""'^, the amplitude of oscillation 
increases considerably. This is just the moment the 
vortex core motion suffers a strong disturb (near hole 
1) and starts to slow down, changing the orbit until 
t ^ 8000J^^, when the core is captured by hole 2. 
Therefore, the kinetic energy of the core decreases 
considerably and a large amount of spin waves is 
released during the instant of the capture process. The 
pinned vortex oscillations also generate more and more 
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FIG. 14: The complete orbit of the vortex core until capture 
by hole 2. The vortex starts its motion in the disk center 
at (20a, 20a). The symbols Xi and X2 indicate the positions 
of antidots 1 and 2 respectively. In the studied process, the 
core follows an almost perfect gyrotropic mode, completing 
five laps around the film before feeling a strong effect of hole 
1, 
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FIG. 15: Mean magnetization < > as a function of time 
and its Fourier transform. 

spin waves due to the rapid rotations of vortex center. 
The amplitude of < > remains large after trap- 
ping but now its oscillations are quite small, meaning 
that the pinned vortex polarization practically vanishes. 

For the sake of completeness, we now discuss a sit- 
uation where the holes are not aligned along the same 
line crossing the disk center. Let us consider a disk 
with radius 20a centered at position (20a, 20a) with two 
small holes {gi = 02 ^9 = 0) displayed at (0, 33a) and 
(12a, 8a). Figure (jl^) shows the orbit of the vortex core 
until its capture by one of the holes. Even in this situa- 
tion (many others, with the holes closer and apart from 
each other were also tested) the main facts remain the 




FIG. 16: (Color online) A picture of the instant that the 
vortex core was captured by hole 2. Note the large amount 
of spin waves generated in this process. Such waves can be 
seen by the out-of-plane fluctuations given by irregularity in 
the spin pattern (green surface) on the fihn. 

same: after being excited by the applied sinusoidal field 
the vortex core starts its gyrotropic motion around the 
disk center. Whenever passing near a hole its dynamics 
changes (similarly to the former case studied, say, both 
holes centered along y-axis, as depicted in Figures 
[T^ , and |l^), and after some revolutions it is eventually 
captured by one of the defects. Analogously to the pre- 
vious configuration with two holes, no core reversal was 
observed in our simulations, neither for smaller (only 1 
spin removed) or larger (8 spins removed) defects. How- 
ever, we cannot strictly state that no other sizes could 
trigger core reversal (with two or more holes), once their 
sizes in our simulations can be only largely varied (keep- 
ing its circular shape), say, while the smallest hole is ob- 
tained by removing a unique spin, the next possible size 
demands we take 4 spin away, and so forth. The main 
point must be recalled: if core reversal could somewhat 
occur with two or more holes, the simplest system where 
this phenomenon would certainly be easiest and clear- 
est observed is that comprising a nanodisk with a single 
non-centered relatively small hole. 

V. CONCLUSIONS AND PROSPECTS 

In this paper we have considered an alternative model 
for circular magnetic thin films containing artificial 
defects and we have also predicted an interesting 
way of reversing the magnetization of the vortex core 
which is completely dependent on its dynamics and 
interaction with a relatively small hole, inserted in 
the sample. Therefore, it should be stressed that the 
present switching process is not stimulated by external 
agents, like usually occur by means of the application 
of magnetic or spin-polarized DC current fields. The 
conditions necessary are summarized by the following 
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FIG. 17: The orbit of the vortex core in a nanodisk with 
radius R = 20a centered at (20a, 20a) containing two holes at 
positions X\ = (0,33a) and Xi = (12a, 8a). As observed for 
several geometries with two defects, the fate of the core is to 
be captured by one of the holes. 

properties: first, only one hole must be present and 
localized at a position near the lateral border of the 
disk; we suggest that the center of the hole must be 
placed at a distance on the order of the vortex core 
size away from the disk border. Second, the hole must 
be small and third, the field required to excite the 
gyrotropic mode should be applied perpendicularly to 
the line joining the disk and hole centers. [Although the 
field is an external agent, we emphasize once more that 
its role is restrict for resonantly excite the gyrotropic 
mode, putting the vortex core in motion. Such a field 
has no direct influence on the reversal mechanism, which 
is completely triggered by vortex- hole interaction]. 
Clearly, the suitable parameters of the applied magnetic 
field, relative hole size, etc, should be probed in actual 
experiments. 

As we observed in our simulations, those special con- 
ditions force the core to move in a circular trajectory 
passing exactly in between the hole and disk borders 
generating a large amount of spin waves, which will also 



develop important role in the switching process (for ex- 
ample, see Fig.(^). Nowadays, such defects-jjuere already 
intentionally introduced in nanostructuresEJO but they 
occupy a relatively large fraction of the material (about 
15-25 percent). Hence, the switching mechanism occur- 
ring during the vortex motion can be met only by fur- 
ther miniaturization of defects in nanostructures. Also ia. 
this line, we should mention the very recent experimentE£l 
where hole defects induce magnetization reversal in elon- 
gated Co rings. The results are aimed at the same di- 
rection of our work and are an experimental indication 
that our proposal may have further and broad relevance. 
An important difference between our predicticp, and the 
phenomenon experimentally observed in Ref.113 is that 
the last needs a permanent presence of an external mag- 
netic field (while in our case the mechanism of the vortex 
switching is based only on internal interactions during 
the vortex motion). In this aspect, the dynamical char- 
acteristic of the phenomenon predicted here is another 
differential with relation to the usual core reversals in- 
duced by external potentials. 

To achieve more details about the vortex-hole interac- 
tions, we have also considered a system with two holes. 
In this case no switching process was observed for several 
configurations analyzed. Since one of the holes will trap 
the vortex, we have then addressed our attention to 
the magnetization dynamics after the capture. When 
the vortex center falls into a hole, some of its energy is 
radiated away in the form of spin waves. The vortex 
center oscillates inside the hole with a frequency much 
larger than the resonance frequency. The bound state 
vortex-hole releases a large quantity of spin waves 
forming coherent magnetization oscillations matched to 
the vortex structure (part of the kinetic energy of the 
core lost in the capture is transformed in spin waves). 
These results demonstrate a significant coupling between 
the vortex (with the core or even captured, without the 
core) and spin waves in a disk. 

To conclude we would like to suggest that the 
switching mechanism reported here can be used in 
technological applications. Actually, this phenomenon 
naturally lends itself to applications in binary data 
storage. Therefore, it seems that our investigations 
may provide fundamentally new ways of using magnetic 
nanostructures in technology. 
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